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RECENT TRANSITS OF THE SUPER-EARTH EXOPLANET GJ 1214B 
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ABSTRACT 

We report recent ground-based photometry of the transiting super-Earth exoplanet GJ 1214b at sev- 
eral wavelengths, including the infrared near 1.25 /mi (J-band). We observed a J-band transit with the 
FLAMINGOS infrared imager and the 2.1-meter telescope on Kitt Peak, and we observed several optical 
transits using a 0.5-meter telescope on Kitt Peak and the 0.36-meter Universidad de Monterrey Observa- 
tory telescope. Our high-precision J-band observations exploit the brightness of the M-dwarf host star at 
this infrared wavelength as compared to the optical, as well as being significantly less affected by stellar 
activity and limb darkening. We fit the J-band transit to obtain an independent determination of the 
planetary and stellar radii. Our radius for the planet (2.61^ is m excellent agreement with the 

discovery value reported by Charbonneau et al. based on optical data. We demonstrate that the plane- 
tary radius is insensitive to degeneracies in the fitting process. We use all of our observations to improve 
the transit ephemeris, finding P = 1.5804043 ± 0.0000005 days, and T = 2454964.94390 ± 0.00006 BJD. 

Subject headings: stars: planetary systems - transits - techniques: photometric 



1. INTRODUCTION 

A current frontier in extrasolar planet research 
is the detection and atmospheric characterization of 
super-Earth exoplanets orbitin g nearby M-dwarf stars 
(|Charbonneau fc Demine 20071). Th e recent detection of 
GJ 1214b (ICharbonneau et alJ l2009h is a significant step 
toward this frontier. GJ 1214b was found by the MEarth 
surve y of nearby M-dwarfs (Nu tzman fc Charbonneaul 
2008j), and it is the first super-Earth exoplanet discov- 
ered by a ground-based photometric method. The ra- 
dius of GJ 1214b derived from the discovery observa- 
tions (2.68ii®, see Table 1) indicates that the planet 
has a low molecular weight composition such as an 
oce an planet, but alternative compositions remain poss i- 
ble (jMiller-Ricci fc Fortnevl[20Tol: iRogers fc Seagerll2010D . 
Inferences concerning the composition of GJ 1214b depend 
on a securely measured radius. In this regard it is valu- 
able to measure the transit in the infrared (IR) . IR transit 
observations are less sensitive to stellar activity and limb- 
darkening, compared with optical wavelengths. Reduced 
sensitivity to these effects helps establish a more robust 
planetary radius. GJ 1214 is also the target of several near- 
term observational programs using both HST and Spitzer, 
so recent ground-based transit observations can benefit the 
space-borne programs by improving the precision of the 
transit ephemeris. 

In this Letter we present results from high-precision pho- 
tometry of a transit of GJ 1214b observed in the J-band 
(1.25 /xm), and several less precise transits in the opti- 



cal. Our IR transit data provi de an independent check on 
the planetary radius derived bv lCharbonneau et a l. (2009) 
from the optical discovery observations, and the totality 
of our recent transit data allow us to significantly improve 
the transit ephemeris. 

2. OBSERVATIONS AND PHOTOMETRY 

2.1. Observations 

We observed transits of GJ1214b using the Universi- 
dad de Monterrey Observatory (UDEM) telescope on UT 
dates April 29, June 6 & 17, 2010. UDEM is a small pri- 
vate college observatory having Minor Planet Center Code 
720, and is located at 689 meters altitude in the suburbs 
of Monterrey, Mexico. The UDEM data were acquired 
using an 7 c -band filter on the 0.36-meter reflector, with a 
1280xl024-pixel CCD camera at 1.0 arc-sec per pixel scale. 
We observed an additional transit on UT May 29, 2010 us- 
ing two telescopes on Kitt Peak. We used the 2.1-meter 
reflecto r with the FL AMINGOS 2048x2048-pixel infrared 
imager (|Elstorjll99§ ). and a J-band (1.25 /mi) filter, at 0.6 
arc-sec per pixel scale. Following the conclusion of nightly 
public programs, we have access to the 0.5-meter telescope 
at the Kitt Peak Visitor Center (VCT) . Simultaneous with 
the 2.1-meter J-band observations, we observed the May 
29 transit using the VCT and a 3072x2048 CCD camera at 
0.45 arcsec per pixel. The VCT observations were planned 
to use the z-band filter, but due to a filter wheel error these 
data were actually acquired in a B filter. Observations at 
all three telescopes used a defocus to improve the pho- 
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tometric precision, and all used off-axis guiding to main- 
tain pointing stability. Exposure times were 30-seconds at 
the 2.1-meter, 60-seconds at the VCT, and 120-seconds at 
UDEM. All of the optical CCD exposures were binned 2x2 
to facilitate rapid readout. 

Flat-field observations were acquired at all three obser- 
vatories using either twilight sky (UDEM and VCT) or 
a series of night-sky FLAMINGOS exposures, incorporat- 
ing pointing offsets to allow removal of stars via a median 
filter. 

2.2. Photometry 

Subsequent to dark current subtraction and division by 
a flat-field frame, we performed aperture photometry on 
the target star and comparison stars. The 20-arcmin field 
of FLAMINGOS provided 8 comparison stars of compa- 
rable IR brightness to G J 1214. We used 7-pixel circular 
aperture radii to measure the stars. We measured the sky 
background surrounding each star using an annulus with 
an inner radius of 9-pixels and an outer radius of 44-pixels. 
Normalizing GJ 1214 to the the comparison stars yielded 
a transit light curve with an observed scatter that varied 
from 0.002 to 0.0015 as a function of time, due to the de- 
creasing airmass during the observations. We found that 
the best results were obtained by averaging the ratios of 
GJ 1214 to each comparison star. This produced smaller 
scatter than the method of ratioing GJ 1214 to the sum of 
all the comparison stars. We estimated an error for each 
GJ 1214 photometric point as the standard deviation of 
the ratio to the individual comparison stars, divided by 
the square root of their number (error of the mean). The 
scatter in the photometry is approximately 65% greater 
than these estimates, so we increased our error estimates 
by this factor to facilitate more accurate % 2 analyses (see 
below) . 

After normalizing to the comparison stars, the GJ 1214 
data clearly reveal the transit but with low-order curvature 
in the out of transit baseline. It is probable that this base- 
line curvature is caused by differences in the effective wave- 
lengths of the J-bandpass as a function of stellar color, in 
combination with the wavelength- variation in telluric wa- 
ter vapor absorption. We account for this baseline using 
a 4th-order polynominal fit, simultaneous with the transit 
fit (see Sec. 3). Photometry of the UDEM and VCT data 
used similar procedures as for the FLAMINGOS data, ex- 
cept that the baseline was modeled as a linear function 
of airmass. Figure 1 shows the 5 transits that we have 
observed. Figure 2 replots the J-band transit, showing 
both the raw photometry and baseline-removed photome- 
try, and overplots a theore tical transit curve based on the 
(jCharbonneau et al.ll2009h system parameters. 

3. DERIVATION OF STELLAR AND PLANETARY RADII 

Our J-band transit data have greater photometric preci- 
sion than our optical data (see Figure 1). Moreover, stellar 
limb darkening and activity signatures are reduced at IR 
wavelengths compared to the optical. Therefore we use 
only our J-band data to fit for the exoplanet and stellar 
radius, but we use all of our transit data to update the 
transit ephemeris. 

Given an estimate of the stellar mass, three unknowns 
(stellar and planetary radii, and impact parameter) follow 



from three constraints (transit depth and duration, and 
ingress/egress time). The solution is only weakly affecte d 
by errors in the adopted stellar mass ([Brown et al.ll200"T[ ). 
Given also that the IR limb darkening is weak, we fix 
both the stellar mass and limb darkening in our fitting for 
radii . We adopt the ste l lar m ass (0.157 solar masses) given 
bv lCharbonneau et al.l (|2009T ). and the model-atmosphere- 
based square-ro ot limb darken ing coefficients for M-dwarf 
stars derived by IClaretj ([1998D . 

We calculate transit curves numerically, using a tile- 
the-star procedure with adaptive-mesh tiles to maximize 
the precision. The trans it curve that results from the 
Cha rbormeau et ail ([2009) system parameters, and the 



Claret ( 19981) limb-darkening, is illustrated in the lower 
panel of Figure 2. For Figure 2, we overlaid this tran- 
sit curve on our data by first removing the baseline from 
the data. We then varied the central transit time to 
find the best fit, not changing the depth or shape of the 
curve. It provides an excellent fit, note particularly that 
it accounts for the sharp ingress/egress portions of the 
data (Fi gure 2). This agreement immediately suggests 
that the iCharbonneau et al.l (|2009T> radii are not signfi- 
cantly affected by limb-darkening errors, or stellar activ- 
ity. Note also that the transit curve still contains signifi- 
cant curvature, unlike tr ansit curves at far-IR wavelengths 
(Ric hardson et al"]|2006t ). where stellar limb darkening is 
virtually zero. The linear limb darkening coefficient for 
this M-dwarf st ar is about a factor of two lower in J-band 
than in I-band (IClaretlll998l). 

Given that the ICharbonneau et al.l (|2009D system pa- 
rameters are very close to best-fit values for our J-band 
data, we identify a refined solution as follows. We 
generate a 2-D grid of transit curves centered on the 
ICharbonneau et "all (|2009f) values, as a function of the stel- 
lar radius and impact parameter. At each grid point, we 
find the best corresponding exoplanet radius by scaling 
the depth of the transit curve and fitting it to the raw 
photometry simultaneously with the baseline polynominal, 
using linear regression. (Scaling the depth of the transit 
curve relies on the excellent approximation that the transit 
depth is primarily dependent on the ratio of planetary-to- 
stellar radius, and not on the impact parameter.) We also 
shift each transit curve in time for the best fit. At each 
grid point we calculate the x 2 of the best fit, and the lo- 
cation of the minimum x gives our global best-fit stellar 
radius, impact parameter, and planetary radius. We de- 
termine errors on these quantities via the range of values 
contained within a y 2 contour of a given significance. The 
results are given in Table 1, except for the best-fit cen- 
tral transit time - that is included with the optical transit 
times in Table 2. 

To fit our four optical transits, we use the lClared(fl998l) 
limb darkening coefficients to generate synthetic transit 
curves at the optical wavelengths. We shift these in time, 
without varying their depth, to find the best central time 
for each transit. Those times are listed in Table 2. We 
combine our Table 2 transit times with th e individual tran- 
sits given by ICharbonneau et aLI ([2009) to derive an up- 
dated ephemeris, listed in Table 1. 
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Figure 3 shows our contours of x 2 f° r the fit to our J- 
band transit, converted to standard deviations, and shown 
in the plane defined by stellar radius and impact param- 
eter. Our x 2 minimum (square on Figure 3) occurs at a 
smaller st ellar radius and sma l ler im pact parameter than 
found by ICharbonneau et al.l ()2009[ ). We do not regard 
this difference as physically significant because: 1) it can 
reflect a degeneracy in the fitting pr ocess, 2) the differ- 
ence is not much greater than the ICharbonneau et al.l 
(2009) error range i n com bination with our errors - the 
ICharbonneau et ail ([2009i) values are within our 2a con- 
tour, and 3) our radius for the planet is in excellent agree- 
ment with ICharbonneau et"aTI (|2009l) . 

The contours in Figure 3 trace a parabolic locus wherein 
a larger stellar radius can trade-off versus a larger impact 
parameter (i.e., transit at higher latitude) to produce the 
same transit duration. In the limit of a small planet and 
zero limb darkening, a chord at high latitude across a large 
star would produce the same transit curve as a chord at 
low latitude across a small star. This mild degeneracy 
is broken by the duration of ingress and egress, so that 
stellar radii exceeding approximately O.23i?0 lie outside 
contours of acceptable \ 2 . We derive a value for the plan- 
etary radius (R p ) by examining the R p values at points 



lying within the la contour (that first contour is not illus- 
trated on Figure 3). This yields R p = 2.61±£?5-R©, where 
the asymmetric error range reflects the fact that the radii 
are not distributed symmetrically within the contour. We 
also derived best-fit values for the stellar radius and impact 
para meter, and these are l i sted i n Table 1 for comparison 
to the ICharbonneau et al.l (|2009T ) results. Our error range 
for the impact parameter (a factor of two) is relatively 
large, reflecting the degeneracy discussed above. However, 
our radius for the planet has a much smaller fractional 
variation within the la contour , and i s in excellent agree- 
ment with ICharbonneau et al.l ()2009f ). We conclude: 1) 
our value for the planetary radius is relatively insensitive 
to the mild degeneracy between stellar radius and impact 
parameter, and 2 ) the discovery value for the radius of 
this super-Earth (Charbo nneau et al.l l2009h is unlikely to 
be significantly affected by uncertainty in the optical limb 
darkening, or affected by stellar activity. 

We are extremely grateful to Dr. Dick Joyce for con- 
siderable and repeated help with FLAMINGOS, that con- 
tributed significantly to the success of our J-band photom- 
etry at the 2-meter telescope. We are also grateful to an 
anonymous referee for a swift and insightful review. 
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Fig. 1. — Photome try of GJ 1214 during 5 tr ansits using three different telescope/wavelength combinations. Nominal transit curves are 
overlaid, based on the Charbonn eau et al.l (120091 ') system parameters, but shifted to fit the best central transit time. 
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Fig. 2.— Top Panel: Raw J-band photometry of the GJ 1214 transit observed from Kitt Peak using FLAMINGOS on May 29 UT. The 
red line is a baseline fit using a 4th order polynominal, necessary due to the different effective wavelengths of the target and comparison 
stars when averaged over the filter bandpass. Bottom Panel: Photometry divided by the baseline, with a theor etical transit c urve (blue line) 
overlaid. The trans i t curv e is based on IR limb darkening for M-dwarf stars (T = & logg = 5.0) from IClareti 1(19981 ). and using the 

Charbonncau ct al. (2009) values for the system parameters. It was shifted in time to give the best fit, but was not altered in depth or shape. 
These J-band data are publically available in the NStED data base. 
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Fig. 3. — Contours of constant \ 2 f° r the best fit to our J-band transit, in the plane defined by the stellar radius and impact parameter. 
The contours have been co nverted from x 2 to s tatis tical significance, and are plotted for 2, 3, 4, and 5cr significance. The cross shows the 
values and ±lcr range from Charbonncau ct al. (2009), and the square indicates the position of our \ 2 minimum. Note the substantial range 
within the contours where the same goodness of fit is obtained by trading a smaller stellar radius for a lower crossing latitude (i.e., smaller 
impact parameter). This is due to a mild degeneracy in the transit fit, as discussed in Sec. 4. 



Transits of GJ 1214b 



7 



Table 1 

Best-fit values for planetary and stellar radi i, impact parameter, and e phemeris parameters, compared to the 

results of cliarbo nneau et al.i (|2009l ') . 



Parameter Best fitted value Charbonneau et al.(2009) values 



Planetary radius 


^■Dl_ .ll-flffi 


2.678 ±0.13i? e 


Stellar radius 


O^Ol-oloos-R© 


0.2110 ± 0.0097i? o 


Impact parameter 


n m +0.181 
U- 1 '°-0.175 


u.ooi_ 082 


Central transit time 


2454964.94390 ±0.00006 (UTC-BJD) 


2454983.9087558 ± 0.0000901 


Orbital period 


1.5804043 ± 0.0000005 days 


1.5803925 ±0.0000117 



Note: ICharbonneau et "ail (|2009l ) list only individual measured transit times, not an averaged epoch. We quote their most 
precise individual transit. 



Table 2 

Times of GJ 1214b transit center, given as UTC-based BJD. 



UT Date 


Telescope 


Filter 


BJD k Error (la) 


April 29, 2010 


UDEM 


Ic 


2455315.79343 ±0.00042 


May 29, 2010 


KP 2-m 


J 


2455345.82126 ± 0.00014 


May 29, 2010 


KP VCT 


B 


2455345.82133 ± 0.00037 


June 6, 2010 


UDEM 


Ic 


2455353.72332 ±0.00036 


June 17, 210 


UDEM 


Ic 


2455364.78669 ± 0.00029 



Note: Adding these transit times to those of ICharbonneau et al.l (|2009f ). gives the best- fit ephemeris parameters (period 
and transit time) listed in Table 1. 



